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1. Introduction 

The quantitative description of soft multiple-parton interaction in hadronic collision has 
been a long-standing and unsolved problem. Such interaction, which for instance con- 
tributes to the underlying hadronic activity, the 'underlying event', has received much 
attention recently, both in terms of experiment at Tevatron Q and theory ||2|, ^, |^, ^, |^. 
With increased collision centre-of-mass energy, one expects a proliferation of such hadronic 
activity. The current operation of Tevatron, as well as the operation of LHC being in plan 
for the near future, adds urgency to its study. 

In terms of experiment, in the simplest language, the underlying event appears as a 
uniform plateau of hadronic transverse energy which must be subtracted from jets. The 
uncertainty in this jet energy correction is the dominant source of the systematic error on 
jet energy at low transverse momentum even at CDF. The effect is expected to become 
more significant ^ with increased centre-of- mass-energy, such as at LHC. 

There has recently been significant progress ^ in the measurement of hadronic activity 
in limited regions of phase space that are defined with respect to the direction of the leading 
jet. These results have been fitted with the available Monte Carlo event generators ||^, |8| 
with some success. 

Another notable measurement has been that of the rate of double parton scattering in 
the process pp ^ j/tt^ + 3j + X at CDF Q . Their measurement gave a surprisingly small 
value for the 'effective cross section' a^s = 14.5 it 1.7^*^2 3 ™b which, however, seems not to 
have been taken account in many theoretical studies since then. 
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Progress has been made on the theoretical side mainly by considering multiple (2 
2) scatters whose frequency is governed by a given 'overlap function'. Events are then 
generated by coupling the (2 — > 2) parton-level scatters with some preferred Monte Carlo 
event generator. 

The superficial convergence of the theoretical effort, which might give one the impres- 
sion that all that is left is to tune the Monte Carlo event generators, is misleading. The 
following are a few theoretical reasons as to why this is the case. 

First, the fragmentation of low-transverse-momentum (Iow-pt) partons is not well- 
understood. Monte Carlo event generators are tuned to describe data from, for instance, 
LEP [10|. Here the usual prescription of perturbative parton shower up to a cut-off followed 
by hadronization, and hence the Monte Carlo event generators, work very well. This is 
because the distribution of hadrons is governed primarily by the partons emitted in the 
parton shower phase and this phase is well-understood. What happens in the low-pj- case 
where there is potentially little perturbative emission is far from clear. 

Second, the above 'overlap function' approach is based on the assumption that the 
proton is a 'disk' with a radius-dependent density, and partons within the proton interact 
independently of each other. This assumption needs justification. In particular, when the 
interaction distance becomes comparable or larger than the proton 'radius', this approach 
would seem insufficient. Even for a process containing a high pT subprocess, the lowest 
energy component of the process occurs near the QCD scale, which is comparable with the 
proton radius. 

Third, the mean number of scatters is proportional to the inclusive cross section which 
is critically dependent on the pT cut-off p™™. At the moment, this imposition of p™™ is an 
ad-hoc procedure and we need to obtain more understanding of the behaviour of the low 
Pt scatters, both in terms of the inclusive cross section and the fragmentation properties. 

In addition, there are some phenomenological insufficiencies which are possibly related 
to the above theoretical problems. We merely outline the problems here. A more detailed 
discussion will be provided later on. 

First, the mean number of scatters is inversely proportional to the above-mentioned 
'effective cross section'. The overlap function approach, if based on the low-energy proton 
form factor, predicts effective cross sections that are much larger than the measured value. 

Second, the predicted underlying event is too soft, i.e., the hadrons carry too little 
transverse momentum. 

Third, the predicted underlying event is too uniform. 

The purpose of this paper is to propose a step towards the solution to the above 
problems within the framework of the multiple-scattering picture. 

We first adopt the measured value of ciefr and set near the Regge-perturbative 
'transition' scale so that, to a good approximation, each scatter can be calculated using 
perturbation theory. The number of scatters is generated according to the Poisson dis- 
tribution. We find that the level of underlying event activity is in agreement with the 
experimental findings. 

We then study the fragmentation scheme dependence using the HERWIG Monte-Carlo 
event generator M. 
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In a recent study |T^, we suggested that the cluster splitting process in the cluster 
fragmentation model ||l2| of HERWIG may be reducible to a modified low-energy as that 
describes 'unresolved' emission below the parton shower cut-off. In accord with the expec- 
tation from this study, the shape of the measured underlying event is described well by 
adopting the modified hadronization scheme and a larger cluster mass cut-off. 

This paper is organized as follows. In sec. ||, we discuss the current status of theoretical 
and experimental study, and outline our interpretation of the discrepancies. In sec. ^, we 
present and discuss the result of our simulation using HERWIG. The conclusions, and the 
outlook regards future application, follow in sec. 



2. Nature of interaction 

2.1 General remarks 

Multiple-parton interaction is often classified as a unitarization effect. The total inclusive 
jet cross section is defined experimentally by counting the number of jet pairs rather than 
the number of events, and theoretically by the ununitarized perturbative calculation for 
a given minimum transverse momentum p™™. For sufficiently large centre-of-mass energy 
^JscM and small p™™, the inclusive jet cross section, either calculated or measured, in 
general exceeds the relevant total cross section. The scattering probability needs to be 
unitarized so that per event, the mean number of scatters n of a 'type' of interaction that 
has inclusive cross section o"™'^ and total cross section 0"*°* is given by: 

inc 

^=^- (2-1) 

Let us now consider, according to some suitable definition, the whole of non-diffractive 
(ND) events. In an ND event we then have on average Ha = ^''a'^/'^nd scatters of the type 
A. 

2.2 The Poisson approach 

The distribution of the number of scatters is not uniquely determined. The simplest ap- 
proach is to assume that the scatters are, neglecting the possible restriction coming from 
the conservation of energy which turns out not to be a significant effect, independent of each 
other. We then have a Poisson distribution of the number of scatters. The normalization 
is such that the cross section for events containing n > 1 scatters is given by: 

cTn = ^e-Vj?S, (2.2) 
n! 

I' being a constant. The inclusive and total non-diffractive cross sections are then related 
by: 

a'^^ = Y,nan, (2.3) 

n>l 

1 — e~'^ 

0"ND = 2-/'^" = 1, ^ND) (2-4) 

n>l 
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so that: 

n = ^ = ^^. (2.5) 

n — > 1/ in the hmit of large i^. Now let us define the following quantity, applicable to this 
case, which matches with the 'effective cross section' of ref. Q: 

tot 

^ _ "ND ^ _tot /n f.\ 

1 — e 

In ref. |^], was measured for the 1.8 TeV run at CDF by comparing the rate of 
back-to-back jet+7/7r'' events that also contain a back-to-back dijet pair, against similar 
events due to pile-up. Systematic uncertainties were hence reduced compared with previous 
studies and they obtain: 

(Teff = 14.5 ± l.rtl'l mb. (2.7) 



We note that this number is small compared with the measured value ||1^ of the nondiffrac- 
tive total cross section at CDF. At 1.8 TeV they find: 

total total total 

^NSD — '^inelastic ~ ^SD 

= (60.33 ± 1.40) - (9.46 ± 0.44) = 50.87 ± 1.5 mb. (2.8) 

Here SD stands for 'single diffractive' events where one of the protons remains intact. NSD 
stands for non-SD. 

The surprisingly large difference between the two cross sections aes and o"^gp can be 
accommodated within the above naive Poisson framework if the remaining 35 mb of the 
NSD cross section is in fact due to the exchange of colour-singlet objects, the Pomeron 
and Reggeons, and does not resolve the quark structure of the proton. If we adopt this 
viewpoint, the 'non-diffractive' cross section should perhaps be understood rather as a 
'resolved' cross section in analogy with the resolved photon cross section. Later on in this 
paper, we use the notation <Jres to indicate this cross section. 

2.3 The overlap-function approach 

A commonly adopted approach in parametrizing the underlying events is to modify the 
picture of independent scatters by introducing an 'overlap function'. In this picture, an 
event occurs at a definite impact parameter b, and scatters are independent per event, i.e., 
eqn. (p. 2]) is replaced by: 



''"^^^e-^Wrf. (2. 



nl 



The inclusive and total cross sections are now given by: 



_1I1C 



i^ND — 

n>l 



i-oo 

Syrian = / v{b)-Kdb^, (2.10) 



^tot 



n>l 



^a„= / l-e-'^(^) rf. (2.11) 
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After this, one normally adopts a factorized form for namely: 



v{h)=a^^A{b). (2.12) 
A{h) is the overlap function which, from eqn. ( 2.1C| ) satisfies the normalization condition: 



A{b)T:db'^ 



1. 



(2.13) 



We can evaluate a^s from A{h). For two distinguishable scatters A and B with small cross 
sections, i.e., (ta,o"_b being much less than ctnD' have: 



0"eff 



A^{b)TTdb'^ 



(2.14) 



Here we are led into a dilemma. If we are to use a physically motivated overlap function 
based on the proton form factor, the resulting cross section is too large. For instance, in 
refs. IQ, ^ we have: 

(2.15) 



A{b) = ^i^^bfKsifib), 

where /i^ = 0.71 GeV^ and Ki{x) are the modified Bessel functions. Substituting this into 
eqn. ( |2.14| ), we obtain 287r/^^ ~ 48 mb, which is too large. 

We add that in terms of theory, the classical picture of parton-filled protons interacting 
according to the parton density given by some form factor does not have strong justification. 
In particular, the typical soft QCD interaction distance is of order of, or even larger than, 
the proton radius. Even events with the hard subprocess at scales much greater than the 
QCD scale in general contain a usually factorizable soft part. 

Let us consider how the simple Poisson distribution of eqn. ( ^.2]) differs from the case 
with a varying overlap function, eqn. (|2.9| ). A simple test case is that of the Gaussian 
overlap function: 



U{b) 



vr 



(2.16) 



A is a dimension-2 constant. With this choice of I'ib), eqn. ( |2.2[ ) can be evaluated analyti- 
cally, and we obtain: 



vr 

nX 



n-1 



^ ml 



m=0 



nX 



^ J_(Aa-7vr)^ 



(2.17) 
(2.18) 



The inclusive cross section satisfies: 



(2.19) 



n=l 
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whereas the total cross section is given by: 



tot 

a 



n=l 



poo 


1 - e--^^') 


'o 





Hence the mean number of scatters per event is: 

n=— + — — (2.21) 

^tot ^ ' 

This last approximation is accurate up to about (n) = 10. 

The largest terms in the series expansion of the exponential are found at the order 
that is close to the expansion coefficient. We therefore see that for small n, eqn. ( 2.18| ) has 



an approximate 1/n behaviour. At large n, this is replaced by a Poisson-like behaviour. 
(Tn decreases monotonically with n, in contrast with the Poisson case where CTn is largest 
at n = n. 

Thus the difference between the two approaches is in principle significant. 

In any case, the deficiencies of the overlap function approach highlights the need to 
rethink the low-energy behaviour of the scattering cross sections. 

The mean number of scatters per event is not too large at Tevatron for reasonable 
p™™ and the measured value of a^s- For instance, in ref. Q, the mean number of scatters 
per event is 2.4. Hence the effect of adopting different distributions is expected not to be 
drastic. On the other hand, for larger centre-of-mass energy, for instance at LHC, such 
effect may become significant. This should be a topic for future experimental studies. 

In addition, the measurement of Cefj on distinct experimental platforms, such as 77 
and 7p colliders, would be an interesting possibility. 

2.4 The CDF 'charged jet' analysis 

In the analysis of ref. |]l|, hadronic activity in the underlying event is studied with respect 
to the leading 'charged jet', which is defined using a non-standard cone algorithm, counting 
only the charged tracks with pT greater than 0.5 GeV and pseudorapidity r/ between — 1 
and -|-1. Jet px is defined as the scalar sum of the pT of the charged tracks within the jet. 

The charged jet construction algorithm proceeds as follows. First, charged tracks are 
ordered in decreasing pT- One then starts with the highest pT track, and define this as the 
seed to form a jet. Going down the list of the remaining charged tracks, tracks which are 
less than Ai2 < 0.7 away from this seed are combined with it. Here Ai? = \/ Ar/^ -|- IS.cfP' 
as usual, with (j) being the azimuthal angle. The new seed is defined to be at the centroid 
of the two objects. This is defined in the rj — cj) space, using pT as the weight. The pT 
of the new object is the scalar sum of the px of the two constituents. This process is 
repeated until the list is exhausted, after which one starts again using the highest px track 
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LEADING JET DIRECTION 




Figure 1: 'Toward', 'away' and 'transverse' regions defined witli respect to tire azimutlial direction 
of tlie leading cliarged jet. 



that remains. After all cliarged tracks have been assigned to jets, with the possibility that 
some jets may contain just one track, the jet construction algorithm is terminated, and the 
leading charged jet is defined as the jet with the largest pT- 

As shown in fig. |l|, three areas in the r] — cp space are defined with respect to the 
azimuthal direction of the leading charged jet. rj is restricted to be between —1 and +1, 
and the three regions are defined according to the azimuthal angle. The 'toward' direction 
is centred around the leading jet direction and is within ir/S radians from the jet direction. 
The 'transverse' direction is between tt/3 and 2tt/3 radians away from the jet direction. 
The 'away' direction is more than 27r/3 radians away from the jet direction. The three 
regions have equal area of 47r/3 in the ij — (p space. 

The quantities considered in ref. |||] , that characterize the underlying hadronic activity, 
are the average scalar sum of pT of charged particles with pT greater than 0.5 GeV, 
and the charged particle multiplicity A^^^hrgi in the three regions described above. These 
quantities were plotted against the pT of the leading charged jet. 

Out of the three regions, the region that is most sensitive to the underlying event is 
the transverse direction. 

2.5 Current status of HERWIG based simulation 

Fig. ^, taken^ from ref. , shows the current status of the fits to the Tevatron data using 
HERWIG-based simulations. 

The default HERWIG curve is very much below the experimental numbers, but this is 
subject to the tuning of parameters. The numbers shown in ref. [||] are in better agreement, 
although still being lower than the experimental result. 

""^We thank the authors of ref. [H for their kind permission to reproduce their figure here. 
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Figure 2: Comparison of 'transverse' hadronic activity, in terms of A^chrg (left) and (right), 
in various HERWIG-based simulations, taken from ref. HERWIG+'Eikonal' model (solid line), 
HERWIG+'Underlying Event' model (solid dashed), HERWIG+'Multiparton Hard' model (dot- 
ted). Experimental data is shown with error-bars. 

The 'multi-parton hard' model is the result of simulation using JIMMY Q. The 
simulation is based on the overlap-function approach. The plot shown has p™™ = 3 GeV. 
The level of hadronic activity is not sufficiently high, and better agreement with data is 
obtained when the inclusive cross section is boosted by lowering p™"^ to 2 GeV. 

The 'eikonal' model is an extension of the JIMMY approach where the gg — > gg 
scatters below = 3 GeV are also generated by modelling the inclusive cross section as 
an exponential function of p^. The gluon structure function is assumed to have a oc l/x 
behaviour. The dependence on p™™ is, according to their claim |^], weakened, but is 
nevertheless present. The numbers drop considerably when p™™ is lowered to 2 GeV. In 
our opinion, this is simply because the inclusive cross section in reality behaves more as a 
power, ~ Pj^^, than an exponential. 

The 'eikonal' model fits the data best out of the three approaches. This is because 
there is not sufficient hadronic activity in the other two approaches and this model adds a 
large amount of 'soft' activity from below p™™. 

As far as the shape of the fit is concerned, there are two remaining problems that are 
not rectified, which are: 

1. When p™™ is adjusted to fit A'^chrg) there is insufficient p?^'^. In other words, the 
charged particles in the underlying event carry too little pT- 

2. The slope at low pT is too steep, i.e., there is too much activity when pT of the 
leading jet is small. In other words, the underlying event is too uniform. 

The above two points seem to indicate two things. 
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"Transverse" PT Distribution (charged) 



CDF 

PT(charged jet#1) > 2 GeV/c .^.^ uncorrected 

theory corrected 

1.8TeV|ii|<1 



Herwig Total 




1 2 3 

PT(charged) GeV/c 



Figure 3: Transverse momentum distribution of cliarged particles {pr > 0.5 GeV,|77| < 1) in 
the 'transverse' region, for leading jet pT > 2 GeV. Experimental data is compared against the 
HERWIG prediction. The dashed curve shows the underlying event contribution. Figure taken 
from ref. Q. 



First, the description of fragmentation in low pT scatters is not adequate. 

Second, the excessive uniformness of the simulated underlying event suggests that the 
idea of employing the cross section below p™™ might actually not reflect the true nature of 
the underlying event. This same point applies to the model intrinsic to HERWIG, which 
also seems to have too much activity implicitly assigned to low pT physics. 

There is a finding in ref. |^] that, in the PYTHIA framework, one needs a double- 
Gaussian overlap function with a 'hard core' in order to fit the data. We note that this 
is possibly another phenomenological indication of the non-uniformity of the observed un- 
derlying event. 

One possibility for large p™™ is that it should be kept high at the scale of the Regge- 
perturbative 'transition', or where colour singlet interactions become more dominant. This 
would seem to be in agreement with the observation made earlier that Ueff is much smaller 
in reality than is supposed in ref. |^]. If is smaller, larger activity can be obtained for 
smaller a™'^. 

In ref. ||l|, in addition to the measurement of A^^^hrg and p?^^, there is also the measure- 
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Figure 4: The low pr colour-octet exchange cross section (left) can be represented as the imaginary 
part of the bubble insertion in the Pomeron exchange diagram (right). 



ment of the px of individual charged tracks in the 'transverse' region. This is reproduced^ 
in fig. ^ and we see that the model intrinsic in HERWIG produces too many soft tracks. 
In ref. a study was made using the other generators, and the result was that PYTHIA 
gives slightly better description of the px distribution although the slope is still too steep. 

One can easily confirm that the excessive softness of the charged tracks shown in fig. ^ 
is quantitatively consistent with our earlier observation that when p™™ is adjusted to fit 
A'^chrgi there is insufficient p^p"^. 

2.6 ' Pert ur bat ive' nature of the underlying event 

We remarked above that one of the reasons for the discrepancies between theory and 
experiment may be that the scatters take place at larger px than is usually supposed, with 
being determined by the Regge characteristic scale. If so, one would expect that, to 
a good approximation, each scatter may be calculated using perturbation theory. 

To justify this claim in terms of the underlying theory, one needs to argue that the 
scatters from the Regge region do not contribute to the underlying event. This may be 
done in two steps by first showing that the Regge inclusive cross section may be suppressed 
and then arguing that the fragmentation in this region contributes less hadronic activity. 

In the perturbative region, for p™™ reasonably small, 0{> 2) GeV, the dependence 
of the inclusive cross section on p™™ is found to be roughly a^^^ ~ (l/p™™)~^. We have 
found numerically that the hadronic activity in default HERWIG also goes as (l/p™™)^~^. 

On the other hand, if the low px scattering cross section is determined by Regge 
dynamics, one may expect a form of the differential cross section that goes as da^^^/dp'^ ~ 
p'^, with e being a suitable Pomeron intercept. Hence the low px cross section is suppressed 
compared with the perturbative behaviour. We do not yet specify the nature of this 
Pomeron here, but our analysis is inspired by the hard Pomeron picture of refs. |14, 

One possible picture which produces this px dependence is illustrated in fig. ^. We 
represent the 'resolved' cross section by the colour-octet gluon exchange cross section, and 



^We thank the authors of ref. M for their kind permission to reproduce their figure here. 
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this is written as the imaginary part of the bubble insertion in the Pomeron exchange 
diagram. 

We write the hehcity and colour conserving gluon-proton amplitude at zero momentum 
transfer as: 

Aoivx, - 9iPx.)[^^ = ■ (2.22) 

We have defined figp as a dimensionless quantity. It is proportional to the product of the 
dimension-(— 1) Pomeron couplings to the gluon and to the proton, e is the Pomeron 
intercept, sq is the Regge characteristic scale and is about 4 GeV^ |14]. As this is an 



overall factor only, the normalization may as well be traded between (3gp and sq. 

In ref. it is argued that the Pomeron coupling to off-shell partons cannot be 

considered to be point-like and that one should introduce a form factor with scale 0(1 GeV). 
Thus the above formulation is valid only below this scale. Above this scale, the form 
factor should suppress the cross section in such a way to reproduce the perturbative result 
cr™'^ ~ (l/p™™)~'^. This picture is in accord with the recent studies p5[| . 

The whole proton-proton amplitude at zero momentum transfer then follows: 

Aubbie(PA.PA. - Pa.PaJ = 2^Pl{Nl (^-i^ j (^-^ j . (2.23) 

The factor 2 in the front is for the two transverse polarizations of the gluon. The Pomeron 
coupling is thought to conserve helicity. We have added the subscript to ^bubble to em- 
phasize the point that we are not calculating the total elastic and inelastic cross section 
but only the 'resolved' cross section, here being identified with the colour octet transfer 
contribution. The pp ^ X cross section is then calculated by the optical theorem: 

a'^^ipp -^X) = J-ImAubbie(pp ^ PP) ■ (2.24) 



SCM 



t=0 



Neglecting the proton mass, let us define the four-momenta as follows: 

PI = ^(1,0,0,1), 
P2 = ^(1,0,0,-1), 

Pg = iE,pT,-p,). (2.25) 

Pz is positive definite and satisfies Pz > \E\. 

We express the phase space in terms of px and rapidity as: 

1 vr 
d'^Pg = -d^pxdtgdy = -dp^dtgdy. (2.26) 

Here tg = P^- We added the subscript to avoid confusion with the overall momentum 
transfer in the amplitude, which is zero. We define rapidity as: 

2/ = -ln^. (2.27) 

^ 2 pz-E ^ ' 
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The two invariant masses squared in eqn. ( |2.23| ) are given by: 



Ml 



Pg + V^cuiPz + E) = tg + J scui-tg - Pt) 



M. 



29 



Pg + ^scuiPz - E) = tg + e scM{-tg - Pt 



As M^g, > 0, we have the inequahties: 



\y\ < In 



^JscM{-tg - Pt) 



(2.28) 



< Pt < -tg - tl/scM- 



(2.29) 



The gluon virtuahty terms tg in eqn. (|2.2g| ) are small and so let us neglect them. This 
is admissible as the cross section is small-x dominated in the structure function language. 
The y integration yields a logarithm which is approximately ln(scM/so)- The remaining 
integration in eqn. ( 2.23| ), expressed in terms of {da^^/ dtg), is merely: 



dt„ 



c 



11 



167r3sot2 



■SCM 

So 



In 



■SCM 

So 



d|m|.| 



So 



(2.30) 



The proton helicities have been averaged over. We have defined m|. = tg -\- which is 
negative definite from eqn. ( p.29| ). The integration is of course trivial. Neglecting the 
small contribution from the lower limit of integration, we obtain: 



da^ 

dtr, 



167r3s2(2 + e) 



So 



So 



V So 



(2.31) 



Thus the cross section has a {—tgY behaviour and falls to zero at zero gluon virtuality 
rather than increasing indefinitely, {—tg) can be loosely exchanged with and hence we 
have demonstrated the power behaviour of the low pj- cross section claimed earlier. 

It is possible to interpret the above discussion in terms of the hard Pomeron picture 
of ref. [|l^, with e ~ 0.4. This provides a means of fixing the normalization of eqn. (2.31), 
by using the gluon structure function. 

Let us consider the first diagram of fig. ^ but with the upper proton replaced by the 
gluon. We write the gg gg amplitude at zero momentum transfer, averaged over the 
first gluon colour and helicity, as: 



99t 



t=o 



iM^dggiM^). 



(2.32) 



The whole amplitude is then: 



■^{9PXp 9PXp, 



2i(3gpiN^ - 1 

MN^-l) 
167r3 



d\ 



(27r)V 



J^lVlfgagg{Mfg) 

9 

2 J„,7\#2 



dtgdp'TdyM(gdgg{M(g) 



SO ) 

(< 

[ So 



1+e 



1+e 



(2.33) 
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We are interested in the structure function fg{x,Q^) = {da/dx)/a. We define x = {pz + 
E)/^scM and omit tg in eqn. ( 2.2^ ) as before, to obtain: 



2 1^ fta=Q'^ J4. / ^2 \ l+e 



figp{Nl-l) dt 



-nin 



^4-41 . (2.34) 

Hence, neglecting the lower limit of integration as before, we obtain: 

^^^"^'^^^16vr3(l + e)(2 + 6) V^o^j ' ^ ^ ^ 



We may compare this with the form given in ref. |15f| : 

xg{x, Q2) ^ 0.95(Q2)1+^«(1 + QVO-5)"^""'/^^-'''. (2.36) 

The units for are in GeV^. This form is claimed to be valid in the range 5 < < 500 
GeV^. eo is the hard Pomeron intercept which, in our notation, is e ~ 0.4. Because of 
the presence of the extra denominator factor (1 + Q^/O.S)^^^'^''/^, the normalization is 
ambiguous. For instance, if we extend their claimed range of applicability to the limit of 
very low and match the two expressions at = 0, we obtain a cross section which is 
unreasonably large. Let us match the two equations at the lowest end, = 5 GeV^, of 



the range of validity of eqn. (|2.36D . We obtain: 



167r3(l + e)(2 + e; 
At 

V^CM — 1-8 TeV, we then obtain from eqn. 



ygp\ c ) (2.37) 



^«30mbGeV-2(^) . (2.38) 
dtg \ So J 

As stated above, there is ambiguity in the normalization due to our lack of knowledge 
about the structure function at very low Q^. 

2.7 Fragmentation in low pr scatters 

In addition to the Regge suppression of low pT scattering rate, we note that fragmentation 
may also be affected by Regge dynamics. 

One characteristic behaviour is Reggeization. Partons that are emitted in the pertur- 
bative phase are coloured objects with colour partners, but in the Regge phase, they may 
become, in some sense, colour singlet. If this is the case, large rapidity gaps would be 
formed per scatter and hence the contribution to the underlying event would be reduced. 

This effect is difhcult to quantify, but since the leading effect is to suppress the low px 
contribution, this provides another justification for neglecting it. 



In ref. |11], we proposed that the splitting of large clusters in the cluster-hadronization 
model [12 1 of HERWIG can be connected to a modified QCD coupling that governs those 
emissions which are considered unresolved in the parton-shower phase. 'Clusters' are 
colour-singlet quark-antiquark units which result from the parton shower process followed 
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by a forced g ^ qq splitting. In HERWIG, they subsequently decay isotropically into 
hadrons according to the phase space weight. 

The effect of this modification is not particularly large for processes occuring at large 
hard-process scale, but it is potentially significant for soft processes such as the underly- 
ing event. Therefore the study of the underlying event, in particular relatively inclusive 
quantities such as the px distribution of charged particles shown in fig. |3|, could provide a 
testing ground for models of hadronization. 

One possibility that follows from the modified picture of hadronization is that because 
of the reduced phase space for parton emission in soft scatters, the resulting cluster size 
is typically larger than in the high pj- processes. This is discussed in the next section, 
together with the possibility that this leads to different individual hadron yields. 



3. Simulation and discussions 



Fitting the experimental data with the simulation tools involves, in general, not a small 
number of tunable parameters. However, they are reducible to three components. First, the 
'effective cross section' deff. Second, the inclusive cross section controlled by p™™. Third, 
the distribution of the number of scatters controlled by the overlap function. Normally 
(Toff is not an explicitly tunable parameter, so that one adjusts p™™ to adjust the amount 
of hadronic activity. The distribution of the number of scatters then affects how much the 
multiple-scattering contribution fluctuates. In addition, one sometimes makes the cut-off 
at p™™ smooth, and this also affects the shape of distributions. 

We now abandon the overlap function and adopt a naive Poisson distribution. We 
hence have only one parameter, p™™, that can be tuned to fit the amount of hadronic 
activity, and have an environment to study the fragmentation scheme dependence. 

Although p™"^ is a tunable parameter, we should not set it too far from the onset of 
Regge physics. A reasonable value may be in the range 1 — 5 GeV and so let us adopt 3 GeV 
for now. For a^s, we adopt the measured value of 14.5 mb. We carry out the simulation 
simply by generating a number of QCD (2 2) scatters on top of one another, without 
considering the constraint from the conservation of energy. 

The experimental data is uncorrected and the theoretical numbers are corrected by 
removing, on average, 8% of the charged tracks. 

In fig. |5|, we show the result of a simple simulation using HERWIG. The inclusive cross 
section is found to be about 27 mb so that about two scatters take place per event on 
average. With this value of p™™, the normalization is found to match that of experimental 
data. We note that in ref. Q, a^s is about 60 mb but the inclusive cross section is larger 
and hence the average number of scatters per event is 2.4, not too dissimilar to the case 
studied here. 

Two hadronization schemes are used. The first is the default cluster hadronization 



scheme of HERWIG. The second is the scheme proposed in ref. [11| which uses a Gaussian 
approximation for low energy as to describe the splitting of large clusters. 

The normalization is smaller for the modified cluster splitting procedure based on the 
low-energy Gaussian as- This is because when very large clusters are split, the modified 
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Figure 5: Comparison of 'transverse' hadronic activity, in terms of iVchrg (above) and p™™ (below). 
Simulation at p^^'" — 3 GeV and acs — 14.5 mb. The two hadronization schemes are explained in 
the text. Experimental data is shown with error-bars. Simulation sample size is one million events. 
Error-bars are not shown. 



procedure results m lighter daughter clusters and hence smaller multiplicity. We note that, 
for the same reason, the p™™ dependence of the underlying hadronic activity is milder in 
this approach. This is welcome as it reduces the theoretical uncertainty due to p™™. 

The two discrepancies with experimental data mentioned in sec. Q are still present to 
some extent, though the slope for low values of leading-jet pj- is now in better agreement 
with experimental data. The remaining problem is that the tracks are still too soft. 

This last point is appreciated in a more quantitative fashion by considering the pT 
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Pj (charged) /GeV 

Figure 6: Transverse momentum distribution of charged particles {pr > 0.5 GeV,|ry| < 1) in 
the 'transverse' region, for leading jet px > 2 GeV. Experimental data is compared against the 
HERWIG prediction for the two hadronization schemes. 



distribution of charged tracks in the transverse region, as shown in fig. ^ 

There is better agreement between experimental data and simulation than in fig. |3|, 
but some discrepancy still remains. The difference between the two hadronization schemes 
is small. The most visible difference is for the point at pT between 0.5 and 1 GeV, where 
the modified hadronization scheme results in fewer tracks, in better agreement with exper- 
imental data. 

One possibility, as to the reason why the pT distribution of charged tracks is too 
soft even with large p™™, has to do with the cluster splitting process. We proposed in 
ref. 1 11] that the splitting of large clusters is due to the parton emission that is considered 
unresolved in the usual parton-shower phase. If so, and considering that the scatters take 
place mostly at very low px, the typical size of clusters from scatters at low pT is expected 
to be larger due to the lack of phase space for parton emission. 

Let us consider the cascade splitting of a cluster with mass Mi^n- From ref. |11], for a 
cluster splitting process — > 12 by unresolved parton emission at p'^ = we have: 



o'-MM. (31) 



and the number of splittings is given by: 



(#spiit)~/olog^. (3.2) 
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Figure 7: Transverse momentum distribution of charged particles {pr > 0.5 GeV,|ry| < 1) in 
the 'transverse' region, for leading jet px > 2 GeV. Experimental data is compared against the 
HERWIG prediction for the modified hadronization scheme with CLMAX set to 5 GeV. 

Here Iq is Cf/tt times the integral of low energy as, and is found to be about 0.5. Com- 
bining the above two equations and considering cascade decay with final cluster masses 
Affinal, we obtain: 

In high pT processes, clusters that are formed before the cluster-splitting process have 
relatively low mass, with some large clusters having masses extending up to 0(10 GeV). 
In this case, the right hand side of eqn. ( p. 3D is about 1.5 for (Q) ~ 0.75 GeV. Hence we 
obtain (Mgnai) ~ 1-5 GeV. On the other hand, in the large Minit limit, in other words 
#spnt ^ oo, we obtain (Mfinai) ^ 2 GeV. 

One simple way of modifying the typical cluster mass as a means to study this aspect of 
hadronization dynamics is to change the HERWIG parameter CLMAX, which is the maximum 
allowed cluster mass. This is imposed even in the case of the modified algorithm. Clusters 
with mass greater than CLMAX are split. We modify this, somewhat arbitrarily, from the 
default value of 3.35 GeV to 5 GeV. We adopt the low-energy as prescription for splitting 
clusters, as otherwise there is no incentive for modifying CLMAX. 

The result, for the pT distribution of charged tracks, is shown in fig. |^, and is in very 
good agreement with CDF data. 

The leading-jet pT dependence of the transverse hadronic activity, shown in fig. also 
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Figure 8: Comparison of 'transverse' hadronic activity, in terms of iVchrg (above) and p™"^ (below). 
Simulation at with CLMAX set to 5 GeV. Experimental data is shown with error-bars. Simulation 
sample size is ten million events. Error-bars are not shown. 

shows good agreement with data. In addition to the improved pT per charged track ratio, 
there is also improvement to the shape of the low leading-jet pT region. 

The remaining difference may be due either to our incomplete description of hadroniza- 
tion or to the inaccurate distribution of the number of scatters. We have tested lowering 

to 2.5 GeV and found that this results in too many soft tracks in fig. ^ so we believe 
that the wrong choice of p™™ alone is not a likely possibility to explain the imperfect fit. 

If the above speculations concerning the fragmentation in low pT scatters are correct, 
one would expect that the identified particle yields may be different. With enhanced cluster 
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hadronization scheme 








HERWIG default 


3748 


18056 


0.208 ± 0.004 


Gaussian as 


2698 


15032 


0.179 ±0.004 


Large clusters 


4334 


14474 


0.299 ± 0.005 



Table 1: The proton to charged pion ratio, for tracks with jr/j < 1 and px > 0.5 GeV, in the three 
hadronization schemes adopted in this study. Sample size is 10000 scatters, p™'" — 3 GeV. The 
errors quoted are due to the Monte Carlo statistics only. 

mass, there is more phase space for decay into heavier hadrons, for instance baryons, 
and hence we would expect increased baryon-meson ratio. This could be confirmed by 
measuring the proton-charged pion ratio Rp/^- 

A simple simulation with a small statistics of 10000 scatters, sampling particles with 
Pt > 0.5 GeV and |?7| < 1 only, yields the numbers shown in tab. ||. We see that in the 
case with CLMAX set to 5 GeV, there is about 50% increase in the proton yield as compared 
to the case with the default CLMAX of 3.35 GeV. 

We note that even the unmodified case leads to substantially increased yield for protons 
compared with yields in high px collision. For instance, the numbers from LEP at the 
pole gives Rp/^ « 1/17 0. 

We should say that some caution is needed as regards the numbers, as the identified 
particle yield for baryons is a poorly described quantity in HERWIG. Hence although 
the above argument does indicate that the proton yield is substantially enhanced in the 
underlying event, the exact numbers for the yield should not be trusted. 

We propose the measurement of Rp/^^, and other hadron yields, in the 'toward' and 
'transverse' regions as a function of the leading jet pT as a way to understand the fragmen- 
tation properties of low pT scatters. 

4. Conclusions and outlook 
Conclusions 

We studied the underlying event in the multiple-scattering picture at Tevatron. 

We started by noting some deficiencies of the current theoretical description of the 
underlying event. We noted the possibility that the scatters affecting the underlying event 
take place at above the Regge-dominated regime and are therefore to a good approximation 
calculable using perturbation theory. We also noted that the description of fragmentation 
in low Pt scatters needs to be improved. 

We carried out our simulation on HERWIG using a naive Poisson distribution for 
the number of scatters, choosing p™''^ = 3 GeV and cJefr = 14.5 mb, and found that by 
modifying the hadronization algorithm as suggested in our previous study and by adopting 
a larger maximum cluster size parameter CLMAX, a good agreement with data is obtained. 

If our picture is correct, one would expect an enhancement in the baryon yield in the 
underlying event. We proposed that a measurement is made of this effect by measuring the 
proton-to-pion ratio as a function of the leading jet pT- A confirmation of this effect, or 
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otherwise, could shed hght both on the nature of the underlying event and on the nature 
of the dynamics of hadronization. 

Outlook 

A possible future improvement to our study may come from a better understanding of 
the distribution of the number of scatters per event, or in other words the mechanism of 
unitarization of scattering probabilities. We desire further experimental study, in particular 
the measurement of a^s on distinct experimental platforms, for instance in resolved 7p and 
77 collisions, and at different -^/icM- The large ^/scM available at LHC will give rise to 
larger mean number of scatters per event. This will be significant aid to the study of the 
distribution of the number of scatters. 

Additionally to the results presented herein, we have made a small-statistics calculation 
for the underlying event at LHC using the same and The flattening-off of the 

transverse hadronic activity occurs at larger trigger-jet pT- For sufficiently large pT, both 
A^chrg and are found to be about 5 times greater than at Tevatron. This coincides with 
the latest results using PYTHIA |l8|. On the other hand, the inclusive cross section 
at LHC is 220 mb according to HERWIG, which is about 8 times that at Tevatron. The 
greater part of the difference is presumably due to the difference in the allowed rapidity 
range, i.e., there are more forward scatters at LHC than at Tevatron. We have made simple 
phase-space estimations to confirm this point. 

The HERWIG sub-version 6.505 has appeared recently. One of the major ingredients 
in this release is an interface to JIMMY. It is a simple matter to modify this HERWIG- 
JIMMY package to generate scatters according to a Poisson distribution as in our approach. 
This merely involves the substitution of the overlap function by a step function, i.e., A{b'^) = 
©(o'cff — '^b'^)/<^eS- We have made some calculations by adopting this, and the results are 
consistent with those presented in this study. 



Acknowledgments 

The author thanks Mike Seymour and Bryan Webber for discussions, and Jay Dittmann, 
Rick Field, Joey Huston, Mario Martinez-Perez and Ming-Jer Wang for answering questions 
related to the details of the CDF charged-jet study. 



References 



[1] T. AfFolder et al. [CDF Collaboration], Charged jet evolution and the underlying event 



proton anti-proton collisions at 1.8 TeV, \Phys. Rev. D 65 (2002) 092002| ; 
R. D. Field [CDF Collaboration], The underlying event in hard scattering processes, in Proc. 
of the APS/DPF/DPB summer study on the future of particle physics (Snowmass 2001) 
ed. N. Graf, eConf C0106 30 (2001) P501 |eCoTi,/ C010630 (2001) P50l| |hep-ph/0201192[] ; 
Data is available online at http:/ /www. phys.ufl.edu/'rfield/cdf/chgjet/chgjet-intro. html . 



[2] T. Sjostrand and M. van Zijl, A multiple interaction model for the event structure in 



hadron collisions, Phys. Rev. D 36 (1987) 201£ 



-20- 



[3] G. Marchesini and B. R. Webber, Associated transverse energy in hadronic jet production, 
\Phys. Rev. D 38 (1988) 341S . 

J. M. Butterworth, J. R. Forshaw and M. H. Seymour, Multiparton interactions in 



[4] 



[5] 



[6] 



photoproduction at HERA, \Z. Physik C 72 (1996) 637| | |hep-ph/9601371 |. 

I. Borozan and M. H. Seymour, An eikonal model for multiparticle production in hadron 
hadron interactions, \j. High Energy Phys. 0209 (2002) 015| ||hep-ph/0207283 |. 

T. Sjostrand and P. Z. Skands, Multiple interactions and the structure of beam remnants, |j] 



High Energy Phys. 0403 (2004) 053 [ iep-ph/040207£ ] 



[7] G. Marchesini, B. R. Webber, G. Abbiendi, I. G. Knowles, M. H. Seymour and L. Stanco, 
HERWIG: A Monte Carlo event generator for simulating Hadron Emission Reactions With 
Interfering Gluons. Version 5.1 - April 1991, Comput. Phys. Commun. 67 (1992) 465; 
G. Corcella, I.G. Knowles, G. Marchesini, S. Moretti, K. Odagiri, P. Richardson, 
M.H. Seymour and B.R. Webber, HERWIG 6: An event generator for Hadron Emission 



Reactions With Interfering Gluons (including supersymmetric processes), I. High Energi 



Phys. 0101 (2001) OlOj |hep-ph/0011363| ; HERWIG 6.5 release note, |hep-ph/0210213 

T. Sjostrand, P. Eden, C. Friberg, L. Lonnblad, G. Miu, S. Mrenna and E. Norrbin, 
High- energy-physics event generation with PYTHIA 6.1, Comput. Phys. Commun. 135 
(2001) 238| [ ^ep-ph/0010017[ . 



[9] F. Abe et al. [CDF Collaboration], Double parton scattering in pp collisions at = 1.8 
TeV, \Phys. Rev. D 56 (1997) 3811 . 

[10] The latest HERWIG parameter settings, tuned with LEP data, are available on line at 



littp://hepwww. rl.ac.uk/theory/seymour/herwig/HWtune. html 



[11] K. Odagiri, Local charge compensation from colour preconfinement as a key to the dynamics 
of hadronization, \J. High Energy Phys. 0307 (2003) 022| |hep-ph/0307026|]. 



[12] B. R. Webber, A QCD model for jet fragmentation including soft gluon interference, Nucl. 
Phys. B 238 (1984) 492[ 



[13] 



F. Abe et al. [CDF Collaboration], Measurement of the antiproton-proton total cross 
section at ^ = 546 and 1800 GeV, \Phys. Rev. D 50 (1994) 555C. 



[14] A. Donnachie and P. V. Landshoff, pp and pp elastic scattering, Nucl. Phys. B 231 (1984) 



189| ; Dynamics Of Elastic Scattering, \Nucl. Phys. B 267 (1986) 690 



[15] A. Donnachie and P. V. Landshoff, Small x: two pomerons!, |P/tt/s. Lett. B 437 (1998) 408| 
[ hep-ph/980634^ ; 

J. R. Cudell, A. Donnachie and P. V. Landshoff, Perturbative evolution and Regge 
behaviour, ^hys. Lett. B 448 (1999) 28l| [Erratum ibid. 526 (2002) 413] [|hep-ph/9901222f ; 



A. Donnachie and P. V. Landshoff, New data and the hard pomeron, Phys. Lett. B 518 



(2001) 63 [iep-ph/0105088|; Perturbative QCD and Regge theory: closing the circle, Phys. 



Lett. B 533 (2002) 277 [iep-ph/0111427|; The proton's gluon distribution, Phys. Lett. B 



550 (2002) 160 [hep-pli/0204165|. 



[16] A. Donnachie and P. V. Landshoff, Exclusive rho production in deep inelastic scattering, 
Phys. Lett. B 185 (1987) 403| ; Charmed quarks: structure function and nonperturbative 



interactions, Phys. Lett. B 207 (1988) 319 



[17] K. Hagiwara et al. [Particle Data Group Collaboration], Review of particle physics, Phys 



Rev. D 66 (2002) 010001 



- 21 - 



[18] S. Mrenna, private communication. 



